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In-medium Yukawa theory is part of the thermodynamics of the Standard Model of particle 
physics and is one of the main building blocks of most effective field theories of fermionic systems. 
By computing its pressure we investigate the nonperturbative thermodynamics at finite temperature 
and density using the optimized perturbation theory (OPT) framework. Our calculations are valid 
. . . for arbitrary fermion and scalar masses, temperature, chemical potential, and not restricted to 

0^ ' weak coupling. The model is considered in the presence as well as in the absence of condensates. 

' Comparison with nonperturbative results shows that second order perturbation theory (PT) fails 

, in the first case but performs rather well when condensates are absent, even at high-temperature 

regimes. 

^ ■ PACS number(s): ll.lO.Wx, 12.38. Cy, 11.15.Tk 



I. INTRODUCTION 



(N 

' (— I ' I Matter under extreme conditions often requires the use of effective field theories in the description of its thermo- 
O ^1 dynamical properties, independently of the energy scale under consideration. For instance, in the phenomenon of 
I ' spontaneous symmetry breaking and temperature dependence of antiferromagnetic order in high Tc superconductors 
described by a field-theoretic version of the Hubbard model [l| at relatively low energies, one has to incorporate the 
f-H interaction of fermions with the crystal lattice via phonons, in a coarse graining of quantum electrodynamics. Moving 
I— —I up in energy scale, the phase structure of strong interactions and the quark-gluon plasma [1], investigated in high- 
energy heavy- ion collisions experiments ^ and in the observation of compact stars [j], often demands simplifications 
[ of quantum chromodynamics, producing a variety of low-energy effective models 

■ Besides the major role it plays in the mechanisms of spontaneous symmetry breaking and mass generation in 
Qv^ ' the Standard Model of particle physics, the Yukawa interaction stands out as one of the main ingredients in the 

I construction of simplifying effective field theories to study the thermodynamics of systems under extreme conditions, 
T-H ' especially if one imposes renormalizability. Although an effective theory does not require renormalizability to be 
, consistent, a physically-motivated cutoff being usually more than satisfactory in this case, this feature can prove to 
' be useful in the study of scale dependence and running via renormalization group methods @ . 
^\ [ In this paper we investigate the full nonperturbative thermodynamics of the Yukawa theory at finite temperature 

■ and density by computing its pressure using the optimized perturbation theory (OPT) framework. In our evaluations 
we consider contributions up to two loops, which include direct (Hartree like) as well as exchange (Fock like) terms, 
so that the Yukawa thermodynamics can be investigated in the presence of condensates and also in their absence. 
In the first case, besides OPT and ordinary perturbation theory (PT) we shall also perform a mean field (MFT) 
evaluation. This is an important step in establishing the OPT reliability since the reader will see that when exchange 
contributions are neglected OPT exactly reproduces MFT results which can be considered "exact" in this (large N) 
limit ^. We present results which are valid for arbitrary fermion and scalar masses, temperature, chemical potential, 
and coupling. The region of large values of the coupling is particularly interesting and useful, since several effective 
field theory models in particle and nuclear physics exhibit Yukawa coupling constants that are much larger than 
one, as is the case, for instance, in. the linear sigma model frequently used in the description of the chiral phase 
transition [l^, [ll|, [l^, [3 IlBl i and in pion-nucleon models extracted from chiral Lagrangians [l^, [l3| . 

When condensates are present we show that, as expected, PT has a poor performance since it cannot resum direct 
(tadpole) contributions associated with symmetry breaking. On the other hand, OPT improves over MFT by also 
incorporating exchange terms in a nonperturbative fashion. 
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^ The reliability of the OPT framework applied to symmetric and broken phases was studied previously in the context of scalar 
theories (see for instance Ref. (7|). 



2 



In the absence of condensates and in the regime of very small coupling our numerical findings can be verified in 
the limit of vanishing temperature by comparison to exact analytic two-loop perturbative results previously obtained 
by some of us for the equation of state of cold and dense Yukawa theory within the MS scheme . In Ref. Q , the 
two-loop momentum integrals were computed analytically for arbitrary fermion and scalar masses, the final result 
being expressed in terms of well-known special functions, which provides us with a solid and clear reference in this 
limit. 

Furthermore, this comparison also provides another way of testing the idea that perturbation theory at high density 
and zero temperature in the symmetric phase is much better behaved than its converse [6, 18, 19i. [20l [2ll. [2^. which 
has well-known severe infrared problems [23j . Our analysis shows that, surprisingly, second-order PT agrees quite 
well with the OPT nonperturbative results up to values of the coupling of order one in the case of cold and dense as 
well as hot and dense Yukawa theory. 

The framework of OPT [24| . also known as the linear delta expansion, is an example of a variational method that 
implements the resummation of certain classes of Feynman diagrams, incorporating nonperturbative effects in the 
computation of the thermodynamic potential (for related methods, see Refs. [1^). It has been successfully applied 
to the study of many different physical situations, such as mapping the phase diagram of the 2-1-1 dimensional Gross- 
Neveu model [26] , where a previously undetermined "liquid-gas" phase has been located, and determining the critical 
temperature for Bose-Einstein condensation in dilute interacting atomic gases p7| . Some early applications at finite 
temperature can be found in Refs. [2^ . 

The paper is organized as follows. In Section II we present the in-medium Yukawa theory and set up the notation 
of the paper. Section III contains a sketch of the derivation of the perturbative thermodynamic potential at finite 
temperature and chemical potential. In Section IV we apply the OPT machinery to the evaluation of the nonpertur- 
bative thermodynamic potential in the Yukawa theory. The thermodynamics of the Yukawa theory, in the presence 
of condensates, is considered in Section V where the results from PT, OPT and MFT are contrasted. In Section 
VI we consider the symmetric case, where condensates are absent. Comparing the results from OPT and PT, we 
show that the latter performs well also for extreme temperature values up to two-loop order. Section VII contains 
our conclusions. Technical details involved in the calculation of the vacuum contributions and direct terms in the 
two- loop thermodynamic potential are left for appendices ^. 



II. IN-MEDIUM YUKAWA THEORY 



In what follows, we consider a gas of Np flavors of massive spin- 1/2 fermions whose interaction is mediated by a 
massive real scalar field, (/), with an interaction term of the Yukawa type, so that the Lagrangian has the following 
general form: 



(1) 



where 



Nf 

^V- = X! - m) tAc , (2) 

>C0 = ^id^<P){d^cj,)-^mlq^^-\:,q^^-Xq^\ (3) 

^znt = V'aV'a'/'- (4) 

Here, m and are the fermion and boson masses, respectively, assuming all the fermions have the same mass, for 
simplicity. The Yukawa coupling is represented by 5; A3 and A are bosonic self-couplings allowed by renormalizability. 
Here we choose A3 = A = disregarding bosonic self-interactions which will be treated in a future work [29j . 

We work in the imaginary-time Matsubara formalism of finite-temperature field theory, where the time dimension 
is compactified and associated with the inverse temperature /3 = 1/T [2^. In this approach, one has to impose 
periodicity (anti-periodicity) for the bosonic (fermionic) fields in the imaginary time r, in order to satisfy the spin- 
statistics theorem. Therefore, only specific discrete Fourier modes are allowed, and integrals over the zeroth component 



^ The technicalities concerning Matsubara sums, renormahzation, and the analytic evaluation of in-medium momentum integrals in the 
calculation of the medium contributions up to two loops were extensively addressed in the appendix of Ref. Q . 
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of four-momentum are replaced by discrete sums over the Matsubara frequencies, denoted by u}^ = 2rmT for bosons 
and Ld^ — {2n + 1)ttT for fermions, with n integer. Nonzero density effects are included by incorporating the 
constraint of conservation of the fermion number via a shift in the zeroth component of the fermionic four-momentum 
= iuj^ ^ = iuj^ + ji, fi being the chemical potential. 
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FIG. 1; Perturbative thermodynamic potential for the Yukawa theory to two loops. Solid lines represent fermions and dashed 
lines stand for bosons. Here we omit the diagrams containing counterterms. The third (tadpole) contribution can be neglected 
when we consider the symmetric case. 



From the partition function written in terms of the euclidean action for the lagrangian ([T]), Zy{T, fi) 
one derives the perturbative series for the thermodynamic potential Qy = ^{^/PV) ^^Zy (23[ : 



Tr exp(— S'y 



e=i 



:t)0 



(5) 



where V is the volume of the system, Zq is the partition function of the free theory and Sint represents the euclidean 
interaction action. Notice that Wick's theorem implies that only even powers in the above expansion survive, yielding 
a power series in ay = g'^/Air. However, at finite temperature the perturbative expansion also contains odd powers 
of g coming from resummed contributions of the zeroth Matsubara mode for bosons, such as in the case of the 
plasmon contribution [23l |. In the zero-tcmperature limit this is not the case and, even including hard-dense- loop 
corrections, only g^^ terms are modified [20]. Since we restrict our analysis to two loops (the resummed OPT 
calculation also departs from a two- loop perturbative setting), we are not concerned with this issue. Omitting the 
diagrams representing counterterms, the perturbative thermodynamic potential is shown, diagrammatically, in Fig. 
[TJ The first two diagrams correspond to the free gas, while the last two represent interaction terms: the third diagram 
is a contribution of the direct type and the fourth is of the exchange type. Notice that the third diagram, which 
contains tadpoles, belongs to the one-particle-rcducible class, that does not contribute in the absence of condensates. 
It will therefore be neglected in Section VI, where only the symmetric case is considered. 



III. PERTURBATIVE THERMODYNAMIC POTENTIAL AT FINITE ^ AND T 

Considering the thermodynamic potential only up to the first non trivial (two- loop) contribution, which corresponds 
to direct plus exchange terms, one must evaluate the following 0{g^) quantity: 



^Y^\j,Jn [K^-ml] +z^Trln[f -,71]-^ ^Tr 



Sdir(0)-HScxc(P) 



(f -m) 



(6) 



where the trace is to be performed over the Dirac structure, and the 4-momenta are given in terms of the Matsubara 
frequencies and the 3-momenta for fermions, P = = ito'^ + M i p); ^-i^d bosons, K — — iujf , k). We choose 
the metric tensor signature g^'^ = diag(-|-, — , — , — ), and the following short-hand notation for sum-integrals: 



The direct contribution is given by 



1 



and, the exchange self-energy is given by the following sum-integral over fermionic momenta: 



1 
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Computing the traces at /Lt ^ and T ^ (see Ref. [y] for details concerning the exchange term), one can write 
the thermodynamic potential as ^ly = Hy'^ + fly'^'^iT, where the vacuum contribution prior to renormalization 
has the form 



(2^) 



3^P 



2mi 



2m 



(Pp 1 



(2^)3 Ep 



Nf f d^pi d^p2 J-{E+~UJi2) 



where cjk -= (k^ + ml^/^, Ep = (p^ + mY^^ E± = Ep, ± Ep,, cj^ = (|pi - P2 



(2^)3(27r)3 LUi2Ep,Ep, 

(10) 

m|)^/^, and we have defined 



J± 



Pl ■P2 ±^p,^p 
~ ^12 



= 1 - 



Aw? — rr?^ 

1?2 ~2~ 

^ ^12 



(11) 



Notice that Eq. pi?)) contains zero-point energy divergent terms which are T- and \i- independent only within the 
PT approach in which case they can be conveniently absorbed by the usual vacuum subtraction which normalizes the 
pressure so that it vanishes at /j, = 0. In practice this means that, as done in Ref. Q, one does not have to care to 
their explicit evaluation and renormalization. However, within the other approaches (OPT and MFT) considered here 
the PT bare mass, m, is replaced by effective (T- and /i- dependent) masses within the same diagrams. Therefore, 
we must renormalize those contributions appropriately. Within the MS subtraction scheme, the fully renormalized 
vacuum term can be written as (see Appendices El and IB] for details) 



REN/ 



-,dir3EN 



(m, m^) + ri' 



Gxc,REN 



(m, m^) , 



(12) 



where we have defined the following functions of the masses: 
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(-)Cxc,REN/_ _ \ 




+6"M-^(l-6^)aK) 



V2 



An? 



7 + log 



m? 



n 2 



V3 



(13) 
(14) 



An? 



(15) 



Here 7 is the Euler constant, A is the renormalization scale in the MS scheme, and the functions Vi{z) and a (n?^ 
are defined in Appendices [X] and |B] 

Adding the direct term to the free gas plus exchange medium-dependent terms evaluated in Ref. [B] one obtains 



rp4 

2^ 



^dz log[l 



'^dz {log[l + e-f^-'^/^)] + log[l 



-(£;,+m/t)i 



g2j,2 F 



(4^4)^2 
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z dz- 
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(2^4)^2 
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2(27r)4 



(2^) 



z^dzy'^dy / du 



^'E^Ey 



z'^dzx^dx 



du. 



t-E+ - iC+E^ 



nb{uj^)N}{l) , 



(16) 
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where, in order to perform numerical investigations, we have defined the following dimensionless quantities: 
+ m^^lT"^, El^ z'^+ m'^/T'^, E±=E^± E,^, El^ = + + 2xz u^^ + mVT^, and 



J± 

Si 
ai 



a2 



Q!3 



{E,TEyY-uo,y^ ' 

n/(ii;,+Ai/r) + n/(£;, -/i/T), 

+ ^i|T) nf{Ey + /i/T) + nf{E, ^ ^l/T) nf{Ey ~ ^/T) , 
nf{E, + fi/T) nf{Ey- /i/T) + nf{E,- /i/T) nf{Ey + /i/T) , 
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(17) 

(18) 

(19) 
(20) 
(21) 



(22) 
(23) 
(24) 



z2 + y2 _|_ 2zy u^y + m^^/T^, x = k/T, z = _Pi/T, y = P2/T, nb{x) = [exp(a;) - 1] \ and 77/(2;) 



with ujly 
[exp(x) + 1]"^ 

Given the perturbative expressions obtained in this section, we can compute the thermodynamic potential fly in the 
OPT framework, which corresponds to resumming all dressed diagrams of the direct and exchange types, generating 
in practice an effective mass for the fermions. 



IV. NONPERTURBATIVE THERMODYNAMIC POTENTIAL IN THE OPT FRAMEWORK 

Following the standard procedure, the OPT framework 24] can be implemented in the Yukawa theory as 



Copt = ipa - + V*)]'^a+S 



^ gff'^ai^a , (25) 



where a sum over flavors is implied, and 77* = 77(1 — 6). As one can easily notice from the deformed Lagrangian above, 
at (5 = the theory is an exactly solvable theory of massive fermions even in the case when, originally, 777 = 0. In 
this particular case, the arbitrary mass parameter, rj, works as an infrared regulator which proves to be very useful in 
studies related to chiral symmetry breaking. We should remark that, due to our choice, the meson sector disappears 
at ^ = since there is no need for their mediation in this case. Finally, when S = I the original, interacting theory is 
recovered, so that our >Copt interpolates between a free (exactly solvable) fermionic theory and the original one. 

The Feynman rules generated by the interpolating theory are trivially obtained from the original ones: 777 1— > 777 + 77*, 
g I— > Sg, and the bosonic propagator receives a 1/(5 factor. It is important to notice that the interpolation docs not 
change the polynomial structure of the original theory and hence docs not spoil renormalization, as proved in Rcf. 
[3O I. Now, a physical quantity such as ily^"'" is perturbatively evaluated in powers of the dummy parameter, S, which 
is formally treated as small during the intermediate steps of the evaluation. At the end one sets i5 = 1, in a procedure 
which is analogous to the one used within the large- TV approximation, where N is formally treated as a large number 
and finally set to its finite value at the end. The result, however, depends on the arbitrary parameterjj^, which can 
be fixed by requiring that ^y^'^ be evaluated at the point where it is less sensitive to this parameter [Srl (Principle 
of Minimal Sensitivity, PMS). This can be accomplished by requiring 



= . (26) 

^1—ri,S—l 



In general, the optimum value 77 becomes a function of the original couplings via self-consistent equations, generating 
nonperturbative results. 
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Now we can apply the method to the case of the two-loop perturbative fiy given in Eq. ([6]). Using the OPT 
replacements for the Feynman rules and expanding 77* to order S, one obtains the first nontrivial result 



OPT 

Y 



S l-T In [K^ - ml] + iT Tr In - (m + 77)] + 5 iT 

^JK JP JP 

■Edir(0,r;) + Eexc(P,??) 



-5 ^T. Tr 



where 



and 



Sdir(0,r7) = ™j( ^) ^Tr- ^ 



m4'J JQ W-im + j])] 



JQ 



1 



(m + 77)][(P-Q)2 -m^l 



(27) 



(28) 



(29) 



Notice that, when compared to the second-order perturbative result, Eq. (I27p displays an extra contribution given by 
the third term on its right-hand side which represents a one-loop graph with a Srj insertion. Now, using Eq. (j26p one 
arrives at 







1 



1 



5 i >r Tr I r „ 



[p - {m + ri)f 



-6 ij: Tr 



Sdir(0,77)-hSexc(P,?7) 



(30) 



where, to obtain the last term, we have used a redefinition of momenta p 1— > g, g 1— > p at an intermediate step. Setting 
i5 = 1, one has a nontrivial, coupling-dependent self-consistent integral relation for the optimum mass parameter 
involving the self-energy given by 



i 



Tr 



?7-Edh.(0)-Sexc(P,j?) 



= 0. 



(31) 



Since the OPT propagator has an infinite number of 77 insertions, one can see that the optimization will resum 
exchange graphs in a nonperturbative way. It is very interesting to notice that when exchange terms are neglected 
77 = l]dir(0) — —5(0)0, where (0)o represents the scalar condensate, satisfying the MET self-consistent relation for the 
effective mass (obviously, when 5—^0 the OPT results agree with the free gas case) . This type of result is consistent 
with applications of OPT to different types of theories [34I and illustrates the way OPT works. Notice that within 
the Hartree approximation one also adds and subtracts a mass term which is determined self-consistently. The basic 
difference is that the topology of this term is fixed from the start: direct terms in the Hartree approximation, direct 
plus exchange in the Hartree- Eock approximation. Within OPT the effective mass m -I- 77 is arbitrary from the start, 
its optimum form being determined by the topology of the contributions considered in the perturbative evaluation of 
a physical quantity such as the thermodynamic potential considered here. 

In the next two sections we shall study the pressure numerically, using different approximations in several different 
situations. With this aim we set ttt, — O.lA^, = 777/2 in our numerical routines The temperature and chemical 
potential ranges considered cover — > IOtti while the coupling values cover — > tt. 



V. RESULTS IN THE PRESENCE OF A SCALAR CONDENSATE 



In this section we compare the results generated by PT, OPT and MET when the scalar condensate represented 
by the direct (one-particle reducible) terms are considered. This case is also interesting because one can use the well- 
established MET to analyze the results provided by OPT and PT. As already mentioned, when exchange contributions 



^ The value = 0, which is of particular interest in situations motivated by QCD, will also be considered in subsection IVI Al 
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are not considered in OPT, MFT results are exactly reproduced since in this case both theories employ the same 
effective mass M^s = m + I]^J^'^(0) — m + fj. Each approximation considers different two-loop contributions: PT 
takes all graphs shown in Fig. [T]into account, while OPT considers all plus the extra one-loop fermionic graph with 
the Srj insertion. In practice, as shown in the previous section, the MFT result is quickly recovered from the OPT 
one by neglecting the exhange term. Figure [2] shows the pressure as a function of /i at T = 0. As one can see, PT 
predicts very high values for the pressure as /i increases, in disagreement with MFT and OPT results. Since OPT 
agrees exactly with MFT when exchange terms are neglected, one can also see in this figure the effects of resumming 
exchange contributions: they yield slightly higher values of the pressure for increasing /i. Figure [3] shows the same 
situation but at a high temperature. In this case, the OPT-predicted pressure values are smaller than the MFT ones 
as /i increases. 
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FIG. 2: The pressure, P/A^jg-, as a function of /i/A^^ at T = 0. Top: results from PT (continuous line), MFT (dashed line) 
and OPT (dot-dashed line). Bottom: differences, due to exchange terms, between OPT and MFT. 



In Figure[4]we analyze the pressure as a function of the coupling for low (0.5m) and high (5r7i) values of /i and T . As 
expected, all methods agree with the free gas case when g ^ 0. Also, at high T, /i values MFT and OPT results tend 
towards the free gas since one approaches the Stefan-Boltzmann limit, while PT has a completely different behavior 
in this situation. Finally, let us study the behavior of the OPT and MFT effective masses as functions of T and /z, as 
shown in Figure [5] for g — n. Both quantities have a quantitatively as well as qualitatively different behavior at small 
T, fi values but, as seen in the previous figures, this effect does not manifest itself in the pressure, probably being 
compensated by the presence (absence) of exchange terms within the OPT (MFT). As T or ^ increases the qualitative 
behavior of both effective masses becomes the same although there are still quantitative differences. These results 
clearly illustrate how PT is not appropriate to deal with this kind of situation since it cannot resum the condensate 
which arises from the interaction between the scalar field and fermions and which is related to symmetry breaking. 
Our application nicely illustrates the reliability of the OPT results since they exactly agree with MFT at large- 
(when exchange terms are neglected), while allows us to improve over this approximation by resumming exchange 
contributions. 



VI. RESULTS IN THE ABSENCE OF A SCALAR CONDENSATE 



Let us now follow Ref. [fi] and neglect the one-particle-reducible two-loop diagrams containing tadpoles. In this case, 
which is relevant for the situation where symmetry breaking is not present, one expects that PT will perform better 
than in the previous case where the tadpole contributions have been considered, especially in the zero-temperature 
limit. Clearly MFT is not applicable in this situation, and we will restrain our analysis to OPT and PT results. At 
this stage, the reader should be convinced that, as explicitly shown in the previous section, the former method is able 
to generate nonperturbative results which can be used to access the eventual breakdown of PT. 
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FIG. 3: The pressure, P/A^jg-, as a function of fi/Ay^ at T = l.OA^j^. Top: results from PT (continuous line), MFT (dashed 
line) and OPT (dot-dashed line). Bottom: differences, due to exchange terms, between the OPT and MFT. 





FIG. 4: The pressure, P/A|j^, as a function of g for fi — T = O.OSAj;^ (top figure) and fi — T = 0.5Af^ (bottom figure). PT 
corresponds to the continuous line, MFT to the dashed line and OPT to the dot-dashed line. 



A. Cold and dense case with — 



In this subsection, the order-S OPT results will be compared to the second-order perturbative predictions of Ref. 
Q. The case = is of particular interest since in this situation one can write the thermodymanic potential 
= rip"^ -I- ri™°'^(T = 0, ^) in a simple and compact analytic form, in terms of the vacuum part 

lim = 4iVFB^™(m)+ lim J7°^^'^™(m, m^) , (32) 



9 



0.25 \ 
0.20 ^ 

0.10 / , 
0.05 / 

0.00 1 



T/Ams-0 



0.0 0.1 0.2 0.3 0.4 0.5 



0.25 



0.20 

Meff/m 

0.10 



/ \ 



0.05 / 

0.00 : 



0.00 0.05 0.10 0.15 0.20 0.25 0.30 

T/Ams 

FIG. 5: The OPT (dot-dashed hne) and MFT (dashed hne) effective masses (m + fj and m + Edir(O)) in units of m for (? = tt. 
Top: Mcfi/m as a function of M/Ajjjg- for T — 0. Bottom; McS /m as a function of T / for jj, — 0. 
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where p^j: = fi^ — m? and u = fipf — log ^ ■ 

Figure [6] shows the pressure as a function of /i for a large value of the coupling, g = tt. As one can see both methods 
predict very similar results. 
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FIG. 6: Pressure in the absence of tadpoles normalized by A^jg- as a function of the chemical potential ^ normalized by A^j^. 
Dots represent the OPT result and the line stands for PT. The fermion mass is fixed at m = O.lA^jg-, and g = "k. 

To analyze the tiny differences, let us define the quantity AP/Pp — \{Popt ~ Pp)\/Pp where Pp and Popt are, 
respectively, the pressures predicted by PT and OPT. Although the numerical discrepancies appear to be rather 
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small, Figure [7| nicely illustrates that AP/ Pp increases with higher couplings and decreases with higher chemical 
potential values as opposed to the case where the scalar tadpole is present. 




FIG. 7: Difference between the perturbative and the OPT pressures in the absence of tadpoles as a function of the coupling g 
and the fermion chemical potential fi normalized by A^jg. The fermion mass is fixed aX m — O.IA^. 



This behavior can be better understood if one analyzes how the OPT effective mass varies with g and fi. Figure [8] 
shows that the quantity (m + ff)/m deviates from 1 as g increases but, contrary to the case where tadpoles are present 
(see Fig. [5|), approaches 1 as ^ increases. 

Therefore, the nonperturbative OPT results support the PT results of Ref. @, at T = 0, when scalar condensates 
are not considered. Notice that even though we have considered in this subsection only the = case our numerical 
simulations show that the agreement between PT and OPT remains valid for ^ 0. For more details concerning 
the effects of the scalar mass the reader is referred to Ref. ral . 



B. Thermal effects 

So far, we have seen through comparison with MFT and OPT that PT does not give reliable results when condensates 
are present but, as seen in the previous subsection, the situation improves in their absence, at least at T = 0. In 
principle, it is not obvious that PT will furnish reliable results at high temperatures, even when only exchange 
contributions are considered. The aim of this subsection is to analyze this situation. Figure [H] compares the OPT and 
PT results for the pressure as a function of /i in two extreme situations: T = and T = 10m — l.OAj^. In the first 
case both methods agree well, as one should expect from the discussion performed in the previous section, but more 
surprisingly is the high-temperature result where OPT also supports PT. 

Like in the previous T = case, the quantity AP/Pp assumes very small values which increase with high values of 
g but decrease with high values of T and/or /i. The OPT effective mass in units of to, (m + ffj/m, also behaves as 
in the previous case, deviating from 1 as g increases and, contrary to the case where tadpoles are present (see Fig. 
E]), approaching \ as ^ and/or T increase. Therefore, it looks like the absence of direct terms means that m :^ m + fj 
at high T and/or /i, so that PT behaves like OPT. Finally, Figure [TOl suggests that even at relatively high T and /i 




{m+T])/m 



1^ o-«'^^^°^/^/\^s°'' 

FIG. 8: The OPT effective mass m + in units of m as a function of the coupling g and the fermion chemical potential ji 
normalized by when tadpoles are absent. The fermion mass is fixed at m = O.lAj^. 
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OPT and PT deviate from the free gas as the coupling increases, while in the case with condensates (see Fig. U ) 
OPT (and MFT) had a better agreement with the free gas as opposed to PT. Therefore, the correct resummation of 
condensates seems to reduce the effects of interactions at high T and /i. 



VII. CONCLUSIONS 



We have investigated the thermodynamics of the Yukawa model at finite temperature and chemical potential by 
evaluating its thermodynamic potential up to the two-loop level which includes direct (Hartree-like) as well as exchange 
(Fock-like) types of contributions. Three different methods have been considered: the usual perturbation theory (PT) 
with its bare mass, the optimized perturbation theory (OPT) with its effective mass given by the PMS variational 
criterion, and the well-known MFT with its self-consistent effective mass. 

We have considered coupling values ranging from g = to g » 1 in situations where the temperature and 
chemical potential ranged from zero to ten times the highest mass value, which is kept fixed. As discussed in the 
Introduction, the Yukawa model usually emerges in the description of various physical situations, ranging from low- 
energy condensed matter phenomena to extremely energetic QCD matter. In our approach, the characteristic features 
of each physical system will be brought about essentially by the specific values of the coupling g and the three energy 
scales: the fermion and scalar masses m and m^, respectively, and the renormalization scale A^|g, which normalizes 
all the quantities in our plots. As observed previously using the perturbative method 6], variations of the masses can 
significantly affect the thermodynamic potential of the Yukawa theory, yielding extremely different thermodynamical 
pictures. The renormalization scale Aj^ sets the typical energy scale of the system of interest. In hadronic physics, 
for example, it is reasonable to choose Aj^ ~ 1 GeV, the confinement scale. In this paper, we keep the discussion in 
general grounds, normalizing physical quantities by Aj^ and fixing the masses as m = 0.1 Aj^ and = , 0.5 to, 
and concentrating on the effects of Hartree- and Fock-like interactions and nonperturbative corrections. 

First, we have analyzed the pressure with both direct and exchange contributions, with the former being associated 
to the presence of a scalar condensate driven by the interactions with fermions. We have shown that OPT and 
MFT are identical if one does not consider exchange terms, in agreement with many other applications [32], which is 
reassuring since in the limit of direct contributions only the MFT resummation can be considered as "exact" . Moreover, 
nonperturbative effects of exchange contributions are readily incorporated by considering OPT consistently up to two 
loops. In the hght of the nonperturbative approaches OPT and MFT, our results show that, as expected, naive PT 
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FIG. 10: The pressure, P/A^jg- as a function of g for ^/ = T = O.OSAjjfg- (top figure) and at = T = O.SAjjfg- when tadpoles 
are absent. PT results are represented by the continuous lines while the OPT results are represented by the dot-dashed lines. 



is inadequate to deal with this situation since it has no ability to resum tadpoles. 

As a byproduct of this application we could see how the resummation of exchange terms performed by OPT corrects 
the MFT framework, which corresponds to the leading order of a 1 /N type of approximation. Having established the 
reliability of the OPT, we have followed Ref. Q imposing the absence of tadpoles at T = 0. We have then shown that, 
in this case, the OPT results turn out to be very similar to the ones given by PT. Finally, still in the limit in which 
condensates are not present, we have investigated the high-temperature case where one could expect the breakdown 
of PT. However, our results have shown that this is not the case and the numerical differences between the OPT and 
PT results are very small. Our results also suggest that, at high T and /i, the presence of condensates minimize the 
effects due to interactions when these contributions are properly resummed. 

One should notice that the direct application of these results as effective-model predictions in different physical 
contexts is restricted to definite energy regimes, within which the relevant physical degrees of freedom can be translated 
into a Yukawa model. For instance, when discussing the thermodynamics of cold and dense baryonic matter, the 
Yukawa model considered here might be a suitable effective theory only at small values of the chemical potential, 
e.g. at < 350 MeV. To describe the properties of baryonic matter at higher values of /i, a different framework is 
necessary to account for the phenomenon of color superconductivity (for a review, see Ref. (ssj). One alternative 
could be to extend this model by adding boson fields from another, non-singlet, multiplet of the color SU{?i) group. 
Nevertheless, since Yukawa-type interactions are almost ubiquituous in the description of fermionic matter, the analysis 
in the present paper provides an understanding of the interplay between direct and exchange contributions to the 
thermodynamics as well as the role played by nonperturbative effects in the scalar Yukawa sector of any given extended 
model. 
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APPENDIX A: VACUUM THERMODYNAMIC POTENTIAL AND RENORMALIZATION 

In this appendix, we address the details involved in the explicit derivation of the vacuum contributions to the 
two-loop thermodynamic potential of the Yukawa theory. In particular, we concentrate on the calculation and renor- 
malization of the 1-loop bubble diagrams and the exchange diagram in the vacuum. The vacuum and in-medium 
direct contributions (the third diagram in Fig. [1]) , which were not considered in Ref. , are left for the next appendix. 

After computing the Matsubara sums (cf. Ref. 6] and the next appendix), the pieces which are not explicitly 
dependent on T and/or ^ correspond to the vacuum contribution in Eq. (jlOp. 

The first two terms in Eq. pH]) can be expressed in terms of the following UV-divergent function 



d3p 



BiM) ^ -\ J dM mJ^ 



1 



(Al) 



Those divergences are cancelled by a field-independent counterterm in the Lagrangian, commonly known as a vacuum 
expectation value subtraction or a cosmological constant. Within the MS subtraction scheme, the tridimensional 
momentum integral above is renormalized to |,13i] : 



B 



REN 



(M) = 



647r2 



3 / A2 



(A2) 



The two- loop T, /i- independent exchange contribution to the thermodynamic potential fly , the last term in Eq. 
T0| . can also be written in terms of UV-divergent vacuum integrals; 



vac 



9 



4(m2 + P-Q) 



(Q^ - m^) {P^ - m^) (Q-P) 



(A3) 



corresponding to the vacuum exchange diagram ^, as shown in Fig. 1111 



Qexc 
"vac 




FIG. 11: Vacuum contribution of the exchange term written in terms of the associated vacuum diagram. 

The renormalization is then implemented through the usual procedure, with the addition of the appropriate diagrams 
containing counterterms, as represented in Fig. 1121 



VAC 



exc,REN 




CVAC ^^Xs^ 



VAC 
/ + 



VAC 



= Qexc I QfCT I QbCT , y 

— ^ 'vac ~r ^ 'vHC ~r ^ 'vac I ^ 



FIG. 12: Renormalized exchange contribution to the vacuum thermodynamic potential. The crosses indicate counterterm 
vertices. 

The counterterm vertices are defined in Fig. [T3] in terms of wavefunction and mass counterterms. At this order 
within the MS subtraction scheme, these counterterm vertices cancel exactly the 1-loop vacuum self-energy poles, 
yielding [d — A ~ e): 



1 



2(47r)2 



si 



2). 



(4^) 



6^1' = - 



(47r)2 



(2) ^ _12m2 2 
(47r)2 e 



(A4) 
(A5) 



* Throughout the appendices, whenever we refer to vacuum diagrams, wo adopt the Feynman rules from Ref. [sj], with factors {—Np) 
associated with fermion loops excluded from the diagram definition. 
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FIG. 13: Definition of counterterm vertices. 



Therefore, the second and third contributions in Fig. 1121 the vacuum bubble diagrams with counterterm insertions, 
yield respectively: 



n 



fCT 



bCT 



Nf , , f (TP 
*^(-2) / Tr 



^^(-2)^(^g2) 4 



2 J (27r) 

.Nf 
~2 



i 



(A6) 



(A7) 



where we have defined the following divergent integral 



drP__}__ 

(27r)4 P2 - m2 



(A8) 



whose regularization yields 



rR-EG/-^2 



(47r)2 



- + 1 + log I ^ ) + e a(m") + O(e^) 



with 



(A9) 



a{rn?) 



If , AA2 

l + loi 



4 I 6 



(AlO) 



Using the results (|A4p . (|A5[) and (jA9p in (|A7|) and (jA6[) . we arrive at the following final regularized expressions for 
the second and third terms in Fig. [T^l respectively: 



2 ^ ^ .y (-4^^4 



bCT 



6m 



- + -logf4)+2«K) + 0(6)^ , 



e e 

2X2 



2 2 , / A2 , „ , 2 
- + -log - > ' 



+ 2a(m^) + 0(e) 



(All) 
(A12) 



In Fig. [121 the contribution still to be calculated explicitly in a regularized form is r^vac- Using the identity 

n?+P-Q = 2n?- \ml -\[{Q - Pf - m^] +\{P^ - n?) + \[Q^ - m") , (A13) 



we can rewrite ri^aci Eq. (IA3p . as 



(A14) 



15 



in terms of the following integrals: 



-'2 



^3 



d^Pd^Q 1 

(27r)8 _ ^2) (p2 _ „,2) [(g _ p)2 _ ^2 

(27r)8 (g2 - m2) (P2 _ ^2) ' 
d^Pd^Q 1 



Defining 



(27r)8 (g2_^2) (g_p)2_^2 
d4pd4g 



J(a,6) = 



(2^)8 (P2_a)(g2_5) 



(A15) 

(A16) 
(A17) 

(A18) 



we have = J{m ,m ) and /| = J(m , m^). The dimensional regularization of J(a,b) is straightforward, yielding: 



3^A2 



{Any \ An 



r 1 



iab)i-\ 



(A19) 



where d = 4 — e. 

On the other hand, the evaluation of the integral if is extremely involved, essentially due to the absence of 
factorization of terms containing different masses. This calculation was performed in Ref. [35| and the result is 



(-')"^^ (I) { - + 2z) + I [4z log(4z)] 



2z 



(2^)2 



[log(4z)]^ + 2(1 - z)$(z) + 0(e)} , 
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4m^ 
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[r(i + e/2)f 

(l-6/2)(l-6) 

3 
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(A21) 
(A22) 
(A23) 



^z) EE 4z { [2-log(4z)] 2i^i ( 1,1,2; ^ 



d { 3 
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dc 



2^1(1, l,c; z) 



and 2-^1 is the hypergeometric function, defined by: 



2Fi(a,6,c; z) = ^ 



fc=0 



{a)k{h)k 
{c)k kl 



(A24) 



(A25) 



where {a)k = ^r(Q)''* is the Pochhammer symbol. In Ref. [s^ . simplified expressions for ^{z) valid in the regions 
z > 1 or z < 1 were also derived. 

Taking the results in Eq. (|A19p and in Eq. (jA20p into the Eq. (|A14p and expanding around e = 0, one obtains, 
after a long algebra: 
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2 647r4 



poles in e = 



7 + log 




2 













Vl 



Am? 



Am? 



, , A2 

7 + log — 



V2 



Am^ 



0(6) , 



(A26) 



16 



where 



vi{z) = 2(70 - 4/32) - 8z (270 + (^2) + 16/32 + 4(1 - z)^ $(z) + 

+ log(4z) {8z [2(1 ~z) f3i+ 71]} + [log(4z)]^ {Az^ - 6z} 
V2{z) = 2(71 - 4/3i) - 8z(27i + f3i) + 16^i z^ + 8z(3 - 2z) log(4z) , 
V3{z) = -6-24Z + 16 



(A27) 
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Finally, the 2-loop contribution to the vacuum expectation value counterterm, X , is then defined in order to cancel 
the poles in Eqs. jJTT]), (|XT2)) and ((X26)) . 

Therefore, collecting the finite parts of Eqs. (jAll|) . (|A12|) and (IA26|) . we obtain the final expression for the 
renormalized exchange vacuum contribution to the thermodynamic potential: 
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APPENDIX B: VACUUM AND IN-MEDIUM DIRECT CONTRIBUTION TO THE THERMODYNAMIC 

POTENTIAL 

The third term in the diagrammatic expansion of the perturbative thermodynamic potential, shown in Fig [U is 
the direct term. In this appendix, we concentrate on the explicit evaluation of both vacuum and in-medium parts of 
this contribution. As usual, the full renormalized form of the direct term of the thermodynamic potential is obtained 
through the addition of the appropriate counterterms, as shown in Fig. [TH 

FIG. 14: Renormalized direct contribution to the thermodynamic potential, represented by in-medium diagrams. As before, 
the crosses indicate (vacuum) counterterm vertices. 
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FIG. 15: Definition of the renormalized tadpole integral. 



Defining the renormalized tadpole integral, as in Fig. 1151 we can rewrite the renormalized direct contribution to 
the thermodynamic potential as: 



dir.REN 



(iV 2 





1 







;'t-REN 
y-^tadpolc ) 



(Bl) 



Following the MS prescription, the counterterm is defined through the renormalization of the 1-loop fermionic 
self-energy, cancelling exactly the pole in e = of the tadpole integral, yielding: 
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Finally, taking these results into Eq. (|Bip . we obtain the following expression for the renormalized direct term of 
the thermodynamic potential: 
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in terms of the dimensionless quantities E'^ = z^ + m? /T'^ and z — p/T 
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